Elastic Compton scattering from the deuteron and nucleon polarizabilities 
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Cross sections for elastic Compton scattering from the deuteron were measured over the laboratory 
angles # 7 = 35°-150°. Tagged photons in the laboratory energy range _E 7 = 84-105 MeV were 
scattered from liquid deuterium and detected in the large-volume Boston University Nal (BUNI) 
spectrometer. Using the calculations of Levchuk and L'vov, along with the measured differential 
cross sections, the isospin-averaged nucleon polarizabilities in the deuteron were estimated. A best- 
fit value of (a — /3)=2.6±1.8 was determined, constrained by dispersion sum rules. This is markedly 
different from the accepted value for the proton of (57— /3) p =10.0±1.5±0.9. 

13.40.Em, 13.60.Fz, 14.20. Dh, 25.20.Dc 



Elastic photon scattering from deuterium can, in prin- 
ciple, yield basic information on the substructure of the 
deuteron and hence the nucleons themselves. Compton 
scattering from the proton has been used extensively to 
determine the polarizabilities of the proton (see Ref. |l]] 
and references contained therein). The electric (ct p ) and 
magnetic (ft p ) polarizabilities constitute the first-order 
responses of the internal structure of the proton to ex- 
ternally applied electric and magnetic fields. 

The current status of the proton polarizabilities has 
been reported in Ref. M. 



[a-ft) p = 10.0 ±1.5 ±0.9, 
(a + ft) p = 15.2 ±2.6 ±0.2, 



(1) 
(2) 



where the first error is the combined statistical and sys- 
tematic, and the second is due to the model dependence 
of the dispersion-relation extraction method. The units 
of a and ft are 10 -4 fm 3 . There is also a dispersion sum 
rule which relates the sum of the polarizabilities to the 
nucleon photoabsorption cross section. The generally ac- 
cepted result for the Baldin sum rule [|| is 



a - 



ft) p = 14.2 ±0.5, 



(3) 



although recent reevaluations yield 13.69 ± 0.14 j|] and 
14.0 ± (0.3 - 0.5) [Q. Note that the experimental value 
(Eq. |^) is in agreement with the sum rule. The polariz- 
abilities, obtained from Eqs. [l] and ||, are 



a p = 12.1 ±0.8 ±0.5, 
ft v = 2.1 T 0.8 T 0.5. 



(4) 



The status of the neutron polarizabilities is much less 
satisfactory. The majority of measurements of the elec- 
tric polarizability of the neutron have been done by low- 
energy neutron scattering from the Coulomb field of a 



heavy nucleus. There is considerable disagreement be- 
tween the two most recent measurements |5|,|| . Schmied- 
mayer et al. Q] reported a value for the static electric 
polarizability of a n — 12.0 ± 1.5 ± 2.0, where the first 
uncertainty is statistical and the second systematic. The 
difference between the static (a) and the Compton (a) 
polarizability is small for the neutron. These data have 
been reinterpreted by Enik et al. , and they have sug- 
gested that a value of a n = 7—19 was more appropriate. 
In a separate experiment, Koester et al. Q have reported 
a value of a n = ± 5. Clearly, the current experimental 
value of a n has large uncertainties. Once a n is obtained, 
ft n can be determined via the sum-rule relation for the 
neutron 101, 



(a + ft) n = 15.8 ±0.5. 



(5) 



The recent reevaluations of this sum rule yield 14.40 ± 
0.66 @ and 15.2 ± (0.3 - 0.5) @. 

An alternate method of measuring the neutron polar- 
izabilities is through the use of the quasi-free Compton 
scattering reaction <f (7, 7' n)p in which the scattered pho- 
ton is detected in coincidence with the recoil neutron. In 
certain kinematic regions, the proton behaves as a spec- 
tator and the scattering is primarily from the neutron. 
There has been one measurement reported on this reac- 
tion using bremsstrahlung photons with an endpoint of 
130 MeV ||. However, due to poor statistics, the re- 
sulting determination of a n effectively gives only an up- 
per limit. Experiments have recently been completed at 
Mainz ||, LEGS, and SAL, but results have yet to be 
reported. 

A third method to determine the polarizability of the 
neutron is through the elastic Compton scattering reac- 
tion d(j,j)d. The only reported measurement of this 
reaction was conducted at Illinois |lC(| . The low photon 
energy of this experiment resulted in reduced sensitivity 



1 



to the polarizabilities and hence large error bars. Since 
the deuteron amplitude is sensitive to the sum of the 
proton and neutron polarizabilities, to obtain specific in- 
formation on the neutron, it is necessary to subtract the 
proton polarizabilities. More serious concerns are the 
contributions from meson exchange currents and other 
nuclear effects. Therefore, the final results for the neu- 
tron polarizabilities depend on a model calculation as 
well as knowledge of the proton polarizabilities. 

The present measurement was performed at the Sask- 
atchewan Accelerator Laboratory (SAL). The facility 
houses a 300 MeV linear accelerator (LINAC) that in- 
jects electrons into a pulse-stretcher ring (PSR), produc- 
ing a nearly continuous wave (CW) electron beam. The 
LINAC and PSR were used in conjunction with a high- 
resolution, high-rate photon tagger, a cryogenic target 
system, and a large-volume Nal detector. Complete de- 
tails of the experiment can be found in Ref. . 

The rf(7, 7)d cross section was measured using tagged 
photons in the energy range 84-105 MeV with a resolu- 
tion of 0.3-0.4 MeV. An electron beam of 135 MeV and 
^65% duty factor was incident on a 115 /im aluminum 
radiator producing bremsstrahlung photons, which were 
tagged via the standard photon tagging technique using 
the SAL photon tagger The average tagged flux 

was ~ 6 x 10 7 photons/s integrated over the photon en- 
ergy range. The tagging efficiency was measured approx- 
imately every 8 hours during the experiment by using 
a lead glass Cerenkov detector directly in the beam to 
detect photons in coincidence with electrons in the focal 
plane. The tagging efficiency was approximately 53%. 

Photons scattered from the 12.7 cm long, liquid- 
deuterium target were detected in the large-volume 
Boston University Nal (BUNI) gamma-ray spectrometer. 
BUNI is composed of five optically-isolated segments of 
Nal, each 55.9 cm in length: the core (26.7 cm in di- 
ameter) and four quadrants that form an 11.4 cm thick 
annulus around the core. Since the inelastic contribu- 
tion to d(7,7')d begins only 2.2 MeV below the elastic 
peak, it was essential that the photon detector have at 
least 2% resolution at 100 MeV. The excellent resolution 
of BUNI is mainly due to the fact that it effectively con- 
tains 100% of the electromagnetic showers created by the 
incident photons. Scattered photons were detected at lab 
angles of 35°, 60°, 90°, 120°, and 150°. 

A zero-degree (or in-beam) calibration of the detector 
was done once, in the middle of the experiment, in or- 
der to obtain both the lineshape of BUNI and an energy 
calibration for the Nal core. The Nal quadrants were 
calibrated daily with a radioactive source (Th-C). 

Figure |l| depicts the BUNI energy spectrum, after ran- 
doms have been subtracted, summed over the incident 
photon energy range. The contribution from empty- 
target backgrounds has also been subtracted. Each chan- 
nel in the focal plane of the tagger corresponds to a differ- 
ent incident photon energy and hence a different energy 



for the Compton scattered photon. To sum over tag- 
ger channels, the detected photon energy in BUNI was 
shifted to a value corresponding to the maximum inci- 
dent photon energy of 105 MeV. The magnitude of the 
shift was determined by which tagger channel registered 
the photon. 
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FIG. 1. Photon energy spectrum for 9^— 150° showing the 
lineshape fit (shaded area) to the elastic peak. The inelas- 
tic contribution is readily visible. The region of interest is 
indicated by the vertical lines. 



After background subtraction, the elastic peak in the 
BUNI energy spectrum was integrated over a specific re- 
gion of interest (ROI), depending on the angle, to obtain 
the yields. This region was chosen to ensure that no 
inelastic contributions were contained in the integrated 
region. The ROI for the 150° photon energy spectrum is 
shown in Fig. |l| (vertical lines). The inelastic contribu- 
tion is evident in the energy range 80-90 MeV. 

Systematic errors associated with the energy calibra- 
tion were investigated by comparing the chi-square be- 
tween the data and the in-beam lineshape for different 
calibrations. At the 95% confidence level, the uncertain- 
ties in the cross sections were 1-4%, depending on angle. 

The detection efficiency for the scattered photons was 
largely constrained by the small size of the ROI and pho- 
ton absorption before reaching the detector. Since the 
ROI (Fig. [l]) only selected the peak of the elastic scatter- 
ing distribution, a significant number of events from the 
tail of the lineshape were excluded. The ROI efficiency 
was deduced using the zero-degree Compton spectrum, 
shifted to the appropriate energy, and represents the re- 
sponse of monochromatic photons in BUNI. Comparisons 
of measurements and EGS simulations of zero degree and 
scattering geometries [p| confirmed that the lineshape 
of BUNI measured at zero degrees was consistent with 
the measured lineshape in the scattering geometry within 
~3%. The ratio of counts inside the ROI to the total 
counts in the lineshape yielded an efficiency of (62±2)%. 

The absorption efficiency was broken into two parts: 
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that due to absorption of photons in the target and as- 
sociated apparatus, and that due to absorption of pho- 
tons in the materials located between the Nal crystal 
and the target enclosure. The first absorption factor 
was obtained with an EGS simulation p3| . The second 
was found by integrating the entire zero-degree lineshape 
from the in-beam calibration and comparing it to the 
incident photon flux as determined by the photon tag- 
ger. The absorption factor was approximately (89±1)%, 
which gave an overall detection efficiency of (55±2)%, 
relatively independent (<1%) of scattering angle. 

The effect of the high-rate photon flux on the measured 
yield was investigated. Rate effects were found to be <2% 
in all cases Jll| . 

Sources of systematic errors in this experiment in- 
cluded target thickness (2.5%), solid angle (1.6%), de- 
tection efficiency (3.6%), incident photon flux (1%), and 
energy calibration (1-4% depending on angle). Adding 
contributions in quadrature gives a total systematic error 
of 6.4% (35°), 4.9% (60°), 4.8% (90°), 4.9% (120°), and 
5.2% (150°). 

The final center-of-mass system (CMS) differential 
cross sections for the d(~/, j)d reaction as a function of 
CMS scattering angle are displayed in Fig. || and are 
listed in Table |. The data were averaged over the inci- 
dent photon energy range of 84-105 MeV. The error bars 
in the figure are the quadratic sum of the statistical and 
systematic errors. 
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FIG. 2. Differential cross sections for elastic scattering of 
84-105 MeV photons from deuterium. The curves are calcu- 
lations of Levchuk and L'vov with different polarizabilities. 



Levchuk and L'vov have calculated the differen- 

tial cross section for the d(~f, ^)d reaction in the frame- 
work of a diagrammatic approach. The scattering ampli- 
tude is expressed in terms of resonance and seagull am- 
plitudes. The resonance amplitudes correspond to two- 
step processes and include rescattering of the interme- 
diate nucleons. The one- and two-body seagull ampli- 
tudes involve a photon being absorbed and emitted at 



the same moment within the energy scale involved. The 
one-body seagull diagrams include the nucleon polariz- 
abilities. The contribution of the various ingredients to 
the d(-y,~/)d reaction at 94 MeV is illustrated in Fig. ||. 
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FIG. 3. Theoretical calculations of Levchuk and L'vov 
for the d(7,7)d differential cross section at E 7 =94 MeV. 
The contributions, added successively, are the reso- 
nant (dotted), one-body seagull without polarizabilities 
(short-dashed), two-body seagull (long-dashed), resonant 
rescattering (dot-dashed), and polarizabilities (solid). 



The solid line in Figs. || and || is the full calculation 
of Levchuk and L'vov with nominal free values of the 
isospin-averaged polarizabilities, a — 12.0, (3 — 3.0 or 
(pi — 8) — 9.0. By imposing the sum-rule constraints 
(Eqs. H and ||) , the difference of the polarizabilities can be 
taken as a free parameter and fitted to the experimental 
data. The dashed curve in Fig. || corresponds to the 
best-fit value of 



(5-/3) = 2.6 ±1.8, 



(6) 



which was determined by minimizing the chi-square 
(x 2 /Nd.o.f. = 2.8/4) between the calculation and the 
data. Choosing the sum-rule results of Babusci et al. 
gives a <4% increase in (a — (3). The fitted value is sub- 
stantially different from the nominal free value, and is 
driven by the back-angle cross sections. The extracted 
(a — [3) may be interpreted in terms of (1) medium mod- 
ifications to the nucleon polarizabilities (Aa and A/3), 
(2) neutron polarizabilities which are quite different from 
the proton, or (3) missing physics in the calculation used 
to extract the polarizabilities. It is important to note 
that the extracted value of (a — (3) is dependent upon 
the theoretical model used. For a discussion of the model 
dependence of the existing calculations see Refs. [p|,pT[. 

Medium modifications to the free polarizabilities for 
Compton scattering from light nuclei have been postu- 
lated Recent measurements on 4 He [l7|] , 12 C [ pi] , 
and 16 O |l8| have suggested that changes on the order of 
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A/3 = -A5 = 4-8 were required in order to describe the 
data. These differences were driven by the back-angle 
cross sections which tended to be larger than the theo- 
retical predictions, as is also seen in the current d{pj,^)d 
results. Other measurements, performed at Lund, have 
reported no modification of the free polarizabilitics [jl9| . 
Clearly this issue has yet to be resolved. In the current 
measurement, modifications of A/3 ^3 would be required 
to explain the data. However, for a lightly bound sys- 
tem like the deuteron, medium modifications might be 
expected to be small. 

Assuming that the proton polarizabilities are unmod- 
ified in the deuteron, the neutron polarizabilities can be 
extracted from the fitted (a — (3). Using the fitted value 
along with Eq. [| yields 



(a-P) n = -4.8 ±3.9, 

which can be used along with Eq. 
neutron polarizabilities, 



a r , 



5.5 ±2.0, 
10.3 T 2.0. 



(7) 

to determine the 



(8) 



The error bars do not include any model dependence 
introduced by the theoretical calculation and are anti- 
correlated due to application of the sum rule. These 
results are surprising, since the neutron values are not 
expected to be radically different from the proton. 

Evidence for modification of the backward spin polar- 
izability (7^) has been reported 20 1. The suggested mod- 
ification would increase the theoretical back-angle cross 
sections (Fig. ||) and hence increase the extracted (a — /3) 
by ~3. However, a recent measurement of d{~{, 7 / p)n re- 
ports no evidence for such a modification of 7^ ]2l| . 

There has also been a recent calculation of the d(pf, j)d 
differential cross section within the framework of baryon 
chiral perturbation theory [g2[ . The results of the calcu- 
lation of Beane et al., which includes 0(Q 3 ) terms and 
only a partial set of the 0(Q 4 ) corrections, are similar 
in magnitude to that of Levchuk and L'vov. However, 
Beane et al. warn that their prediction at 95 MeV has 
considerable uncertainty due to the slow convergence of 
the series. Understanding what physics may be missing 
in the theoretical calculations will have to wait for the 
full 0(Q 4 ) treatment. 

TABLE I. Center-of-mass differential cross sections for the 
elastic scattering of 84-105 MeV photons from deuterium. 
The first error is statistical and the second is systematic. 



(deg) 
35.0 
60.0 
90.0 
120.0 
150.0 



(d 7 eg) 
36.8 
62.7 
93.0 
122.6 
151.5 



dcr cm 7dn 7 
(nb/sr) 
13.3 ± 1.3 ± 0.9 
13.5 ± 0.8 ± 0.7 

11.5 ± 0.6 ± 0.6 
15.7 ± 1.1 ± 0.8 

18.6 ± 1.1 ± 1.0 



Other theoretical calculations for the d(j,j)d reaction 
have been reported. The calculations of Karakowski and 
Miller |^3| underpredict the present measurement by a 
factor of ~2 at backward angles when using the nominal 
free values of the polarizabilities. Wilbois et al. [Q have 
reported a calculation at 100 MeV but with polarizabil- 
ities set to zero. Finally, Chen et al. p^j have reported 
results only up to a photon energy of 69 MeV. 

Forthcoming results on the quasi-free Compton scat- 
tering from the deuteron, with detection of the recoil- 
ing nucleon, should shed some light on this controversy. 
By doing the experiments in quasi-free kinematics (in 
the energy region E 7 =200-300 MeV and backward an- 
gles for the scattered photons) the model dependence is 
minimized, and the polarizabilities can be extracted sep- 
arately for the proton and the neutron. 
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